Superconducting proximity effects in metals with a repulsive pairing interaction 
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Studies of the superconducting proximity effect in normal conductor/superconductor (N/S) junc- 
tions almost universally assume no effective electron-electron coupling in the N region. While such 
an approximation leads to a simple description of the proximity effect, it is unclear how it could be 
rigorously justified. We reveal a much more complex picture of the proximity effect in N/S bilayers, 
where S is a clean s-wave BCS superconductor and N is a simple metal with a repulsive effective 
electron coupling. We elucidate the proximity effect behavior using a highly accurate method to 
self-consistently solve the Bogoliubov-deGennes equations. We present our results for a wide range 
of values of the interface scattering, the Fermi wave vector mismatch, the temperature, and the ratio 
g of the effective interaction strengths in the N and S region. We find that the repulsive interaction, 
represented by a negative g, strongly alters the signatures of the proximity effect as can be seen in 
the spatial dependence of the Cooper pair amplitude and the pair potential, as well as in the local 
density of states near the interface. 

PACS numbers: 74.45.+C, 74.78.Fk, 74.78. Na 



I. INTRODUCTION 

For nearly eight decades it has been recognized that 
superconducting properties can leak out from a super- 
conductor into a neighboring metallic region 1 - - — which by 
itself would not be superconducting. This phenomenon is 
known as the superconducting proximity effect. That su- 
perconductivity can penetrate into a non-superconductor 
for a long distance^ has fascinated the condensed matter 
community ever since this was discovered. 

The main signatures of the proximity effect include the 
penetration of the Cooper pairs, with associated phase 
coherence, into the non-superconducting region, and the 
suppression of the pair potential (the superconducting or- 
der parameter) in the superconductor, near the interface. 
Important insights in the proximity effect 4,5 are provided 
by its connection to the process of Andreev-Saint James 
reflection^— An incident electron approaching a normal 
metal/superconductor (N/S) interface from the N region 
can be reflected as a hole, resulting in the transfer of a 
Cooper pair into the S region. This is a phase-coherent 
scattering process in which the reflected particle carries 
information about both the phase of the incident parti- 
cle and the macroscopic phase on the superconductor. 12 
Thus, Andreev reflection is responsible for introducing 
phase coherence in the normal region. Since this reflec- 
tion is a two-particle process, it is plausible to conclude 
that the proximity effect will be also weaker whenever 
this anomalous reflection is suppressed, as e.g., in low- 
transparency N/S junctions. 

Very impressive advances in the fabrication of 
superconducting junctions (including atomically-flat 
interfaces^) have in recent years stimulated extensive 
experimental and theoretical studies of the proximity 
effect. For example, many recent efforts have focused 



on elucidating proximity effects in junctions including 
ferromagnets 1 ^— or superconductors with unconven- 
tional (non s-wave) pairing symmetry ^ 19 ' 20 However, 
significant challenges remain even for the studies of the 
proximity effect in a simple N/S junction, where N is 
the normal conductor and S is a conventional BCS su- 
perconductor with phonon-mediated s-wave pairing sym- 
metry. One such issue is that of the role of the effective 
pairing interaction in the N material. In the simplest 
BCS version of the theory, the S region is characterized 
by a coupling constant A conventionally taken as posi- 
tive for the attractive case. In nearly all of the standard 
treatments of the S/N proximity effect this constant is 
assumed to vanish in the N region^ This implies that 
the pair potential, which enters in the underlying mi- 
croscopic equations, would completely vanish in the N 
region (although the pair amplitude would not) for any 
choice of the N region and the N/S interface. Yet, this 
zero coupling assumption is hardly realistic: while the 
low-frequency phonon mediated interaction is, on general 
grounds, always attractive, the coupling A represents the 
difference between this attraction and the Coulomb pseu- 
dopotential, which is invariably repulsive. The balance 
of the two quantities may lead to a positive A, leading to 
superconductivity, or a negative A but it is most unlikely 
that the two would exactly cancel. Indeed one would 
expect, in non-superconductors, negative values of A in 
roughly the same absolute value range of those found in 
superconductors. 

This was already noted a long time age 22 in a review 
article by P. G. de Gennes, and it was implicit in even 
earlier workj 7 ^ The idea was extensively followed up at 
the time: the comprehensive review article^ mentioned 
above discussed and reviewed many important aspects 
of proximity effect phenomena, with emphasis on issues 
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that could be tackled at the time. The contemporary 
constraints were both experimental and theoretical, im- 
plied by the quality of the samples and interfaces then 
available, the absence of suitable high-resolution probes 
such as the scanning tunneling microscope, and the lim- 
ited capacity of the existing computers, which largely 
restricted theory to analytic methods. Among other top- 
ics, in Ref. [2jthe effect on T c of an attractive or repulsive 
interaction in the N material was considered, with more 
emphasis on the attractive case; a version of the so-called 
"Cooper limit"— argument, for dirty superconductors 
and thin N and S layers, already given in Ref. was 
presented; and a qualitative account of the energy gap 
behavior at an N/S interface was included.— As further 
follow up to these reviews, additional work, including e.g. 
several experiments^— that used electromagnetic and 
transport properties to estimate the pairing interaction 
and its sign, were subsequently performed in the early 
seventies. Yet, despite the wide dissemination of these 
reviews and the high reputation of their authors, activ- 
ity on this problem eventually dwindled and relatively 
little attention has been paid for many years to the ques- 
tion of the influence of a negative value of A in the TV 
material, on the proximity effect. The work of Refs. 
and KM restricted to the quasiclassical limit, and that of 
Ref. |3l| on critical currents in the Cooper limit are among 
the few exceptions. The unspoken assumption elsewhere 
seems to have been that provided A is non-positive its 
value does not matter. Yet, as pointed out already in 
Ref. [22l . such an assumption is quite unreasonable. A 
repulsive interaction in the N region will tend to break 
the Cooper pairs coming from the S region, which are 
responsible for the proximity effect. The situation might 
be in some ways reminiscent of that found in the fer- 
romagnet /superconductor (F/S) proximity^ effect for a 
weak ferromagnet, where the range of the effect is re- 
duced and signatures are found in the local density of 
states (LDOS) due to modified Andreev reflection and 
Andreev-Saint James states^ 

In this work we reexamine this long-standing question 
in a rigorous way. We do this by solving the relevant mi- 
croscopic equations, the Bogoliubov-de Gennes 3 ^ equa- 
tions, in a fully self-consistent way, using computational 
methods recently developed and applie d 33 ' 34 to study 
several aspects of the F/S proximity effects in clean sys- 
tems with smooth interfaces. We consider an N/S bilayer 
in which each layer is thicker than the superconducting 
coherence length, £ 7 in the S material, and we study the 
properties of the system, focusing particularly on both 
the Cooper pair amplitude and the pair potential as a 
function of position, as well as on the LDOS near the 
interface. We study the problem for several values of 
the interface scattering, the Fermi wave vector mismatch 
between the two materials and, most important, of the 
value of the nonpositive ratio g between the effective pair- 
ing constant Xn on the N side and the positive value As 
on the S side. We find that the proximity effect markedly 
depends on g, with the penetration of the pair amplitude 



into N being reduced as the absolute value of g increases. 
There is even a "negative" proximity effect: the presence 
of a repulsive interaction in the normal metal depletes the 
pair amplitude in S near the interface. There are also sig- 
natures of the g dependence on the LDOS as measured 
very near the interface. Not surprisingly, these signatures 
are different from those apparent in the quasiclassical 29 
approximation, which cannot be accurate at very small 
length scales^ As a useful byproduct of our computa- 
tions, we find that it is erroneous, in the study of the N/S 
proximity effect, to subsume the separate effects of inter- 
facial scattering and wave vector mismatch into a single 
effective parameter. This has been long known 3 ^— to be 
the case for F/S interfaces but the situation in the N/S 
case was still unclear. 



II. METHODS 

To study this problem, we solve self-consistently the 
Bogoliubov-deGennes (BdG) equations,— the relevant 
microscopic description for a clean system. The geom- 
etry we consider consists of one normal metal slab of 
thickness d s juxtaposed to a similar slab of an ordinary 
BCS superconductor of thickness ds- We assume a flat 
interface of arbitrary transparency, infinite in the x- and 
y-directions, while the z-axis is normal to the interface. 
The methods we use here have been presented and dis- 
cussed elsewher o 33 ' 34 i 38 ~ — and details need not be given 
again here. In this geometry the BdG equations can be 
written as, 

in terms of the spin- up quasi-electron, u\{z), and spin- 
down quasi-hole, v^z), amplitudes. Here A(z) is the pair 
potential (order parameter—) which is to be determined 
self-consistently as explained below. The single-particle 
Hamiltonian is 

H = k 2 z /2m + e ± + U(z) - E F (z), (2.2) 

where k 2 /2m and e± — fc^/2m, denote the kinetic energy 
from motion in the z and x — y direction for parabolic 
bands, respectively, U(z) is a scalar potential and Ep(z) 
represents the Fermi energies (band widths) Efn, Eps 
in the N and S regions. We set K = ks = 1, for Planck's 
and Boltzmann's constants. The variables e± are decou- 
pled from the z direction, but they affect the eigenvalues 
e n . These variables are measured from the chemical po- 
tential. As pointed out already in many places) 33 ' 34 ' 37 " — 
one should not assume in this type of problem that the 
Fermi wave vectors kpjy and kp$ (or, equivalently, Efn, 
Eps, as measured from the bottom of the bands) are the 
same in both materials. Thus, we introduce a dimension- 
less mismatch parameter defined as 

A = E FN /Eps = (k FN /k FS ) 2 . (2.3) 
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We characterize the mismatch by this single parameter 
without the additional introduction of different effective 
masses in the ./V and S regions,— which would not alter 
our findings. In Eq. (|2.2I) we choose the scalar potential 
to describe the interfacial scattering U (z) — HqS(z — zq), 
where zq is the location of the interface. The strength 
of this scattering is conveniently given in terms of the 
dimensionless barrier strength Hb = mHo/kps, which 
in the limit of no mismatch (A = I), coincides with the 
parameterization introduced in Ref. llO . 

The BdG equations (|2.ip above must be solved to- 
gether with the self-consistency condition: 

= + <(*)vU*)] ^nh , 

n 

(2.4) 

where X(z) = \n < and A(z) = A5 > in the N and 
S region, respectively, the prime in the summation sign 
indicates that it is limited by the usual Debye cutoff, and 
T is the absolute temperature. 

The procedures to diagonalize the BdG numerically 
while ensuring full self-consistency^ have been explained 
in the previous work mentioned above. Basically, one 
chooses a convenient set of orthogonal functions, in this 
and most cases it is appropriate to take sine waves, and 
expands Eq. (I2.1j) and (|2.4I) in terms of that set. The 
required matrix elements are the same as those used, 
for example, in Ref. |4l] in the appropriate (nonmag- 
netic) limit. One assumes an initial form for the function 
A(z) and then iterates the process until self-consistency 
is achieved, that is, the A(z) function obtained from 
Eq. (|2.4p is the same as the one input in Eq. (|2.ip at 
the previous step. 

Once self-consistency is achieved, one can directly ex- 
amine quantities such as the pair potential A(z) and the 
Cooper pair amplitude F(z) (also known as the conden- 
sation amplitude) where, 

F(z) = A(z)/A(z). (2.5) 

It is also very useful to examine, as a more accessi- 
ble experimental quantity, the local density of states 
(LDOS) which is obtained from tunneling experiments, 
where the corresponding spectroscopic information can 
provided by scanning tunneling microscopy. We can ex- 
press the LDOS N(z, e) directly from the self-consistently 
calculated amplitudes (u n<Tl v na ) for the BdG equations 
(EH) as, 

a cr.n 

(2.6) 

One can integrate N(z, e) over a wide region of z to obtain 
the global DOS or over a very small region to obtain local 
results. Since our methods are free of any quasiclassical 
assumptions, our results are reliable even when the region 
examined is microscopic, i.e. of the order of the Fermi 
wavelength. 
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FIG. 1: (Color online) The spatial dependence of the nor- 
malized (see text) pair amplitude F(Z). The dimensionless 
coordinate Z is in units of the Fermi wave vector in the S 
region (Z = kpsz). The interface at Zq = 200 separates the 
N region on the left with repulsive interaction (g — —1/3) 
from the S region on the right which has attractive super- 
conducting coupling. The results are given for temperatures, 
expressed in terms of T c , the transition temperature of bulk 
S, of T/T c of 0.01, and 0.1 through 0.8, from top to bot- 
tom in the S region. The interfacial scattering Hb = 1 and 
the mismatch parameter A = 0.5 [Eq. (|2.3|) ] correspond to a 
low-transparency junction. 



III. RESULTS 

The results of our calculations are described in detail 
in this section. We will measure all the lengths in units 
of the Fermi wave vector hps m the S region, and define 
the relative dimensionless coordinate Z = kpsz. The 
thicknesses of the N and the S regions are taken to be, 
in these units, = D$ = 200 while the superconduct- 
ing coherence length (in the same units) is So = 100. 
Since Ds = 2So, the pair potential rises to very close to 
its bulk value at the far edge of the superconductor. On 
the other hand, it does not decay all the way to zero in 
N: recall that the simplistic estimate^ of the proximity 
depth is of order Tp/T (Tp is the Fermi temperature) 
in a clean system. It is obviously impractical for a nu- 
merical calculation to make the N slab thicker than this 
value, and it is not necessary either for our study, which 
focuses largely in the region near the interface. The val- 
ues of g == Xn/^s studied are g = 0, —1/3, —2/3, —4/5. 
The values of the temperature are given in terms of the 
ratio T/T c where T c is the bulk transition temperature 
in S. Most of the data presented here are at relatively 
5 lowjterppjeratures (T = 0.1T C ) but for most of the values 
of g we have obtained also results for reduced tempera- 
tures of 0.01, 0.2, 0.3 and, in a few cases, up to 0.8 at 0.1 
intervals. Values of the mismatch A [Eq. (|2.3[) ] of 0.25, 
0.5, 1, 2, and (in a few cases) 4 have been studied, while 
for the barrier parameter Hb we have considered values 
of 0, 0.5, and 1. 

In Fig. [Tjwe show that the behavior of the pair ampli- 
tude F(Z) as a function of temperature is as expected. 
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In this figure we have used values of Hb = 1 and A = 0.5 
which correspond to strong interfacial scattering and high 
mismatch: hence a reduced proximity effect, as the N 
and S regions are weakly coupled. The value of g is inter- 
mediate, g = —1/3. F(Z) is normalized so that its value 
in bulk S material at T = would be unity. Results for 
temperatures from nearly zero to 0.8 T c at approximately 
equal intervals are shown. The range of dimensionless 
distance from the left edge Z = Kps z includes one co- 
herence length on the N side (at the left) and nearly all 
of the superconductor. The interface is (in all figures) at 
Z = Zq = 200. We can see that F(Z) rises towards the 
appropriate bulk value deep in S. Because of the strong 
scattering and high mismatch, the profile of F(Z) near 
the interface is rather abrupt and the proximity effect 
overall, quite week compared with other cases discussed 
below. We see also that the temperature dependence of 
the proximity effect is in this case not very drastic [except 
as to the overall level of F(Z)] but appreciable. 
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scription (Hb, A) of the N/S interface, we show in Fig.[3l 
results for the effect of the mismatch parameter A at 
nonzero g (g = —1/3 in this case) at Hb=0- This time 
the quantity shown is F(Z) (at smaller \g\ the N side the 
curves for A(Z) are harder to see). Results are shown for 
five values of A ranging from 1/4 to 4. It is evident that 
the proximity effect is dramatically enhanced when the 
absence of interfacial scattering is combined with the ab- 
sence of mismatch (A = 1), i.e., in the Cooper limit. 
Indeed, in this case (the (red) continuous curve) the be- 
havior of F(Z) at the interface is very smooth, without 
any sign of of abruptness. In the scale shown in the fig- 
ure, and even at the finite temperature (T — 0.3T C ) in 
the figure, the pair potential appears to settle into a con- 
stant on the N side. This is of course not true: it keeps 
decaying but very slowly since the proximity depth is, in 
this case 2 - much longer than the range shown in the figure 
and, indeed, longer than the numerical sample size. 
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FIG. 2: (Color online) The spatial dependence of the pair 
potential A(Z) (normalized in the same way as F(Z)) for 
different values of the interfacial scattering Hb- From top to 
bottom, the results correspond to decreasing values of Hb = 1 
(blue), 0.5 (green), and (red). A(Z) is calculated at low- 
temperature, T = 0.1T C , in the absence of mismatch, A = 1, 
and for a strong repulsive interaction in the N region, g = 
-2/3. 

In the next figure, Fig.[2l we discuss the influence of the 
interfacial scattering parameter Hb at a relatively high 
value of | g \ (g — —2/3) and, for clarity, in the absence 
of mismatch (A =1). The quantity plotted this time is 
the pair potential A(Z) (the "order parameter") as a 
function of Z. A(Z) is normalized in the same way as 
F(Z). Results for the three values of Hb studied are 
shown. In contrast to F(Z), the pair potential in the N 
region will be negative for a repulsive interaction (g < 0) . 
The negativity of A(Z) in N causes it to abruptly jump 
at the interface. We see that this jump increases with 
Hb- This is as one would expect since higher barrier 
scattering isolates the S from the N material and leads 
to increased pair potential in S near the interface and to 
less leakage of Cooper pairs on the N side. 

In a complementary way, for this two-parameter de- 



FIG. 3: (Color online) The spatial dependence of the pair 
amplitude F(Z) for five different mismatch values A =1/4 
(purple), 1/2 (green), 4 (cyan), 2 (blue), and 1 (red), from 
top to bottom on the right side. The results are given for 
H B =0, g= -1/3 at T — 0.3T C . 

In the previous two figures, Figs. [2] and 02 we have con- 
sidered separately the influence of the parameters Hb 
and A. This is actually necessary: it is often assumed 
that these two parameters can be subsumed into a sin- 
gle parameter Z e g that characterizes the effective barrier 
strength. In our notation, Z e g^ would be related to Hb 
and A as, 



Z e fi = 



Ik 

A 1 / 2 



(1 _ A V2)2 



4AV2 



1/2 



(3.1) 



It is instructive to relate this Z e g to the normal state 
junction transparency, i.e., the transmission coefficient 
of an N/N junction at normal incidenc o 10 ' 43 ' 45 



T NN = 1/{1 + Z 2 cS ), 



(3.2) 



implying that Z e g = (or, equivalently, Hb = and 
A = 1) correspond to a completely transparent N/N 
junction, while Z Q g » 1 would correspond to a very 



5 



low-transparency (tunneling) limit. While such a single 
parameter (Z e g) interface description offers a simplified 
approach and has been widely used, it may lead to even 
qualitatively incorrect trends, as compared to the correct 
description in which the effects of Hb and A are consid- 
ered separately. This has been discussed in the context 
of F/S junction o 37 i 46 where it was possible to find an ex- 
ample in which the mismatch can actually enhance the 
junction transparency 36 ^ 3 ' 47 " — 

Our results reveal that Z e g alone is also inadequate 
to describe the proximity effects even in N / S junctions. 
Equivalently, just characterizing such junctions with ei- 
ther the corresponding transmission or reflection coeffi- 
cient, as frequently done^, does not provide an unam- 
biguous description. This can be clearly seen in Fig. [4] 
where two chosen pairs of (iJs,A) values [(0,f/2) and 
(0,2)] both yielding the same Z e s = 0.174 (an interme- 
diate value also equivalent to Hb = 0.174 and no mis- 
match, A = 1) produce strikingly different results for 
both the spatial dependence of the pair amplitude and 
the local DOS. The two curves for F{Z) display very dif- 
ferent suppression near the interface in the S region and 
different decay in N. Analogous differences could also be 
observed for other pairs leading to the same or similar 
Z c ff. Clearly, the difference of calculated pair amplitudes 
which correspond to the same Z c g also implies that other 
quantities such as the pair potential and the LDOS must 
also be inequivalent for the same Z e g. This nonequiva- 
lence is clearly shown for the LDOS near the interface in 
the second panel of Fig. [4] for the same pairs of values 
(0,1/2) and (0,2). There we show the LDOS averaged 
over a region five Z units wide centered at Z = 205, i.e., 
near the interface. The energy is in units of the bulk gap 
in S, Ao, and the LDOS is normalized to the value that 
it would have in an equivalent region in bulk S mate- 
rial: that is, the quantity plotted would be constant and 
unity in the normal state of bulk S material. These nor- 
malizations will be used in all LDOS plots below. The 
inequivalence is obvious. The results have very different 
peak positions and their heights are not at all the same. 
We find that this is true in general. Only in a very crude 
sense are mismatch increases (that is, values of A differ- 
ent from unity) equivalent to increases in Hb, in that 
both lead to diminished proximity effect. With a little 
reflection, and perhaps some hindsight, one can realize 
that, since the nature of the surface scattering (normal or 
Andrccv) originating from the mismatch is not the same 
as that arising from the barrier, the failure of this naive 
approximation should not have been so unexpected. 

Returning now to the basic question of interest, the 
effect of the strength of the repulsive interaction in the 
N region, we show in Fig. [5] results for the pair ampli- 
tude F(Z) for different values of the strength parameter 
g. For the results in this figure the values of Hb = 
and A = 1 correspond to a transparent N/N barrier [re- 
call Eqs. (|3.ip and (13.21) ]. i.e. a very strong proximity 
effect. Another signature of this high transparency junc- 
tion is the lack of Friedel oscillations, seen in the low- 
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FIG. 4: (Color online) Comparison of results for two dif- 
ferent pairs of interfacial scattering and mismatch values, 
(Hb, A) = (0,l/2) and (0,2), which lead to the same effective 
barrier strength Z c a [Eq. (|3.1|) ]. The first panel shows A(Z). 
The top (red) curve at large Z is for the first set of values, 
the other (green) curve is for the second set. The right panel 
shows the local density of states (LDOS) N(E) near (see text) 
Z = 205. The energy E is in the units of the bulk zero- 
temperature superconducting gap Ao. N(E) is normalized to 
its value in the normal state of the bulk S material (see text). 
The curve with the higher (red) peaks is the for first set of 
values, the lower (green) peaks are for the second set. The 
results in both panels, evaluated at T — 0.1T C and g = —1/3, 
clearly show that Z e ff alone can not describe the proximity 
effect or LDOS in N/S junctions. 



transparency case in Fig.[T] F(Z) does not vanish even in 
the farthest region shown on the N (left) side and indeed, 
as expected, shows no sign of decay in the length scales 
shown. Towards the extreme right of the figure, nearly 
two coherence lengths in S, F(Z) approaches in all cases 
its bulk value. However, looking in the regions near the 
interface and in the N region itself, we can clearly see 
how the proximity effect is strongly affected by g: the 
pair amplitude decreases with increasing \g\ not only in 
the entire N region, but also in the S region within over 
one coherence length from the interface. Thus, as stated 
in the Introduction, the influence of a negative g per- 
vades not only the N material but also a thick region 
in the superconductor, where the Cooper pair density is 
rather severely depleted. 
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FIG. 5: (Color online) The spatial dependence of the pair 
amplitude F(Z) for different values of the repulsive interac- 
tion in the iV region: g = 0, —1/3, —2/3, —4/5, from top to 
bottom. The results are given for T = 0.1T C , Hb = and 
A = 1. 

We next examine the LDOS in the middle of the N 
and 5* region, about one coherence length away from 
each side of the N/S interface. The results are aver- 
aged over a region of width 100 centered in the middle 
of the N and S regions (Z = 100 and Z = 300, respec- 
tively). Energy and LDOS are normalized as explained 
in connection with Fig.|4j Panels (a) and (b) in Figure [6] 
show the LDOS evolution with g, for the same parame- 
ters as used in Fig. [5j In both the N (panel (a)) and S 
(panel (b)) regions, N(E) changes smoothly with g with- 
out the appearance of any new features, as compared to 
the g = limit (in the absence of any repulsive interac- 
tion). In this high transparency limit (Hb = 0, A = 1), 
there are strong proximity effects as can be seen in F(Z) 
from Fig. [5j We therefore expect and find a (/-dependent 
LDOS even one coherence length away from the interface, 
in both the N and S regions. We note that the LDOS at 
E = is appreciably larger in the N region. Increasing 
\g\ leads to larger LDOS values near E = 0, while the 
peak near E ± 0.5 moves to slightly larger values of \E\. 
At E = 1 there is a vestige of a peak on the S side, but 
not on the N side. This is of course reasonable, since 
such a peak exists in the bulk S, but not in the bulk 
N, material. In the panels (c) and (d) of Fig. [5] we con- 
sider the effect of the intcrfacial barrier strength on the 
LDOS, averaged in the same way, for a fixed g = —2/3. 
With increasing Hb there is a suppression of the N(E) 
peak at E ss ±0.5 in both the N and S electrodes, but 
an increase in the peak at E — 1 in the S side, as the 
proximity effect decreases. Near E — the LDOS is 
enhanced only in the N region. An increase in Hb di- 
minishes the penetration of A(Z) in the N region and 
reduces its depletion in the S region (see Fig. [5]) and this 
leads to the LDOS one coherence length away from the 
interface looking rather more bulk-like as Hb increases. 
This effect is more marked if, in addition to increasing 
Hb, one sets A ^ 1. 

The effect of the repulsive interaction should be most 



pronounced close to the N/S interface. We focus next, 
therefore, on the spatial dependence of F(Z) and A(Z) 
near the interface (Z = Zo = 200) as a function of g, 
in the low-temperature limit, T = 0.1T C . Panels (a) and 
(b) in Fig. [3 correspond to the high-transparency limit 
with Hb = and A = 1 [recall that Z c s — 0, from 
Eq. (|3.ip j. As seen already in Fig. [3] (the (a) panel in 
Fig. |7| is a blow up of the interfacial region of Fig. [SJ, 
for every value of g F(Z) gradually increases with Z and 
is smooth near the interface. In the limit of g = this 
behavior is well-studied j 21 i 41 On the other hand, right 
at the interface, there is a strong suppression of A(Z), 
as compared to the bulk value. In the S region A(Z) 
and F(Z) decrease very markedly with \g\ reflecting that 
the repulsive interaction in N induces a negative prox- 
imity effect in S. In the N region A(Z) is finite and 
negative for g < 0: this sign change is reminiscent of 
the pair potential behavior due to the formation of 7r- 
states at interfaces with unconventional superconductors 
or F/S junctions^i ^ a ^^ a - H In this region, |A(Z)| 
increases with \g\ and decays away from the N/S inter- 
face. For g < —2/3, the length scale for this decay ex- 
ceeds the superconducting coherence length. In the other 
two panels, (c) and (d), of Fig. [3 we consider, for com- 
parison, parameter values in the low-transparency limit 
(Hb = 1, A = 1/2). In contrast to the other case, there is 
a much sharper, nearly discontinuous rise of F(Z) with a 
step-like behavior near the interface, for every g. The 
rapid (~ kp 1 ) oscillations that can be seen in the S 
region for both F(Z) and A(Z) are not numerical ar- 
tifacts but represent Friedel-like oscillations induced by 
the sharpness of the interface, which would not be seen in 
the quasiclassical treatment of this problem, since such 
an approach would average over a kp 1 length scale. As 
compared to the high-transparency limit, the two main 
differences that can be seen for |A(Z)| in Fig.[3are that 
it attains a much higher value in both the N and S re- 
gions next to the interface; and that it decays much faster 
in the N region, away from the interface. This decay of 
A(Z)| is in general nearly perfectly monotonic in both 
the high- and low-transparency limits. However, we have 
found some cases in which F(Z) has a slight dip just in- 
side the N region, next to the interface. This can occur 
when Hb is large (reflecting an interface scattering po- 
tential, averaged over a Fermi wavelength, of order of the 
Fermi energy) or in the presence of large mismatch (see 
e.g. Fig. [3]). In some cases this translates in A(Z) having 
a minimum slightly away from the interface. An example 
can be seen for A = 2 in Fig. 21 

We explore further the interplay of strong proximity 
effects and the repulsive interaction in the N region by 
considering the interfacial LDOS at a region only one Z 
unit wide centered at Z = 201, just one hundredth of 
a coherence length from the N / S interface, in the S re- 
gion. The results shown in Fig. [8] are given for the high 
transparency Hb and A parameter values used in Figs. [5] 
and [3 Results at Z = 199, just within the N region, are 
very similar. Comparing with the corresponding results 
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FIG. 6: (Color online) The local density of states (LDOS) N(E) [Eq. ijH]) ] in the middle region of the N and the S regions 
(centered at Z=100 and 300 respectively and averaged over 50 Z units (see text) at T = 0.1T C . Panels (a) and (b) show N(E) 
for the N and 5 sides, respectively, at the same parameters used in Fig. [5] Hb — 0, A = 1 and g — 0, —1/3, —2/3, —4/5, 
corresponding to curves (red, green, blue and purple) for which N(0) increases with \g\. Panels (c) and (d) similarly show the 
LDOS for the parameter values used in Fig. [2] g = —2/3, A = 1, and three values of Hb, Hb = 0, 0.5, 1 (red, green and blue), 
LDOS peaks near E — 0.5 are lower with increasing Hb- 



in Fig. |51 we see that the results near the interface are 
qualitatively more similar to those well within (one co- 
herence length) the N side than to those well within the 
S side, and exhibit the same trends. Quantitatively, how- 
ever, there are large differences: the peaks near E = 0.5 
for example, are much more prominent at the interface. 
This is not really surprising: the pair amplitude is al- 
ready drastically suppressed at the interface and (in this 
high transparency case) it remains rather high in N even 
one coherence length away. The increased LDOS near the 
interface arises at least in part from additional Andrccv- 
Saint James states. Interfacial 7r-states have been shown 
to yield low- £7 LDOS peaks as observed in d-wave super- 
conductors and attributed to the formation of Andreev 
bound states (ABS). Furthermore, even for N/s-wave-S 
junctions with a repulsive interaction in the N region, 
low--E LDOS peaks have been predicted using analyti- 
cal but not self-consistent results (a step-function profile 
of A(Z)) or employing quasiclassical approximations to 



calculate A(Z) i 29 i 30 For example, for a spatial profile of 
A(z) qualitatively similar to the one we calculated in 
Fig. O the quasiclassical result of Ref. (the inset of 
their Fig. 1) shows a very sharp zero energy LDOS peak. 
Our results show that, instead, the states at low energy 
are enhanced but the zero-energy peak is absent. This is 
often found theoreticall y 34 ' 41 and experimentall y 59 ' 60 in 
F/S junctions. These results could reflect that the qua- 
siclassical approximation can alter the exact position of 
the low-is 1 interfacial LDOS peak for g < 0£L We will fur- 
ther discuss the Z dependence of the LDOS in connection 
with Fig. Hm below. 

To study in more detail the influence of g on the LDOS 
enhancement in the low energy region that we have ob- 
served, we subtract the corresponding g = LDOS near 
the interface from its g 7^ value. These normalized dif- 
ferences SN(E) are shown in Fig. [51 which corresponds 
to LDOS results averaged over five Z units centered at 
Z = 195. One can readily see that the effects of a 
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FIG. 7: (Color online) The pair amplitude F(Z) and the order parameter A(Z) plotted in regions near the interface for four 
values of g, at T = 0.1T C . Panels (a) and (b) correspond to Hb — 0,A = 1 (high transparency) and panels (c) and (d) to 
Hb = 1,A = 0.5 (low transparency). In all cases lower values correspond to higher \g\. (red, green, blue and purple curves 
correspond to g = 0, —1/3, —2/3, —4/5, respectively) 



nonzero g are quite significant in the low energy region. 
The first panel corresponds to a high transparency junc- 
tion, while the second panel is for a low-transparency 
case. We see that the increase of the LDOS with \g\ is 
largest in the region near E = 0. The maxima are at 
E = in the low transparency case and at nonzero E 
at high transparency. The latter situation has been also 
found to occur in LDOS or conductance results and at- 
tributed to unconventional pairing with broken time re- 
versal symmetry ! 6 ' 37 ' 62 " ^ These peaks for the LDOS dif- 
ference are clearly reminiscent of the peaks for the LDOS 
itself reported in quasiclassical studies for g < N/S 
junctions! 29 ' 30 In the quasiclassical approximation, the 
peak position was studied as function of a reflection coef- 
ficient R that assumedly characterizes overall the specific 
N/S interface. The pitfalls of using a single parameter 
for such purposes have been exposed above, and the use 
of quasiclassical methods to study very narrow spatial re- 
gions at the Fermi wavelength level is obviously suspect. 
It was expected that increasing R would shift the LDOS 
peak position from E — to finite iJ-valuesi^ 9 - Compari- 



son with our rigorous, two parameter, results from Fig- El 
shows the opposite trend for the differential LDOS peak, 
which moves to E = with decreasing transparency. Our 
low-transparency differential LDOS results resemble the 
frequently observed LDOS zero bias conductance peak, 
a signature of ABS in unconventional superconductors, 
which typically becomes narrower with the increase in 
Z c s i.e., in The trend we find, that a low-E the 
interfacial SN(E) increases with g, is rather robust (we 
can see it for both high- and low-transparency and for a 
wide range of temperatures, not just T < 0.1T C ) and it is 
possible that it could be used to directly identify g < 
experimentally. Such an identification would be com- 
plicated by using the suitable g = LDOS background 
subtraction. 

It is helpful to consider in more detail the behavior 
of the LDOS as one approaches the interface. While 
the existence of in-gap states is an inescapable conse- 
quence of the Andreev-Saint James surface states, nei- 
ther experiment^ 9 - nor theor y 33 ] 34 require that they be 
located at zero E. We now examine here how the po- 



FIG. 8: (Color online) The local density of states very near 
the interface (a region one kp S unit wide centered at Z = 201, 
one unit from the interface). At Z — 199 the results are very 
similar. Results in this figure are for Hb = 0, A = 1 and 
(from bottom to top near E = 0) are for increasing values 
of \g\ (red, green, blue and purple curves correspond to g — 
0, -1/3, -2/3, -4/5 respectively). 



sition of the \ow-E LDOS peak depends on the LDOS 
location in the N/S structure and how this dependence 
changes as g varies, and how this might correlate with 
the minimum in A(Z) being pushed away from the in- 
terface, as mentioned above. To elucidate this question, 
we consider the spatial evolution of the LDOS for both 
the g = and g < cases, as shown in Fig. [TDJ Results 
are given for the LDOS, evaluated at regions centered 
at distances 1, 3, 5, and 50 (in the usual units of hp 1 ), 
from each side of the N/S interface, and averaged over 
a region of the same total width as the corresponding 
distance from the interface. As always, we normalize the 
LDOS to the value it would have in an equivalent region 
of the bulk S material in its normal state. In the absence 
of a repulsive interaction (g = 0, panels (a), (b) for the N 
and S sides respectively), moving away from the interface 
reduces the the height of the low- £7 LDOS peaks (below 
E w ±0.5) in both the A^ and S regions. Near E = 
there is virtually no change in the A" region (consistent 
with A(Z) = 0, for Z < 200) while, on the other hand, 
there is a marked decrease in the zevo-E LDOS as one 
moves away from the interface in the S region. 

Comparison of these findings with those at g = —2/3 
(panels (c) for the A^ side and (d) for the S side), re- 
veals a similar LDOS structure but a different situation. 
While there are no new zero-E peaks, there are some 
clear differences. Near E = in the A" region, moving 
closer to the interface does not lead to new LDOS peak 
formation, but we can notice a clear enhancement in the 
zero-E LDOS as the interface is approached from either 
side, while the height of the \ow-E peaks increases very 
markedly as the interface is approached. On the S side 
of the interface, the low E peaks also increase as the in- 
terface is approached, and move slightly towards lower 
energies. We have observed similar trends in the LDOS 
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FIG. 9: Difference between the normalized LDOS, SN(E), 
near the interface (Z — 195) at g = and g < 0. The first 
and second panels correspond to a high- and low-transparency 
interface with Hb — 0, A = 1, and Hb = 1, A = 0.5, respec- 
tively. 



for other junction parameters at other distances from the 
interface. In all those cases studied, we do not find a zero- 
E peak as was associated with ABS in the quasiclassical 
studies.- 9 However, as can be seen from comparison of 
the A~ regions for g — and g = —2/3, the repulsive 
interaction leads to an enhanced \ow-E spectral weight. 
Subtraction of the LDOS, calculated close to the inter- 
face, for 3 = from that of g = —2/3 ((c) and (d) panels 
in Fig. [10]) would again lead to a peak in the differen- 
tial LDOS near E — 0, as shown in Fig. [SJ Thus, the 
observable effect of g is confirmed. 



IV. CONCLUSIONS 

While the possible role of a repulsive effective electron- 
electron interaction in the normal metal of a conduc- 
tor/superconductor (N/S) junction, had already been 
noted in the early seminal work on the superconducting 
proximity effect,— subsequent studies almost universally 
considered such interaction to vanish identically. Per- 
haps the reason was that for the N / S proximity effect 
such a neglect leads to considerable simplifications. In 
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FIG. 10: (Color online) The local density of states for various distances from the interface in the high transparency limit 
(Hb = and A = 1) at T = 0.1T C , Panels (a) and (b) show results for the N and S regions of the sample respectively, at g — 0, 
while (c) and (d) show results for g = —2/3. In all cases the ordering of curves, from top to bottom at E = 0.5 corresponds to 
Z-regions at distances, in units of fc^ 1 , 1, 3, 5, and 50 (purple, blue, green and red curves), from each side of the N/S interface, 
averaged over a region of the same total width as the distance from the interface. 



the N region, the pair potential A(r) (the superconduct- 
ing order parameter) vanishes identically and only the 
pair amplitude F(r) needs to be considered. This leakage 
of Cooper pairs in the N region could then be approxi- 
mately inferred by simply considering Andreev reflection 
and a step-function pair potential, — although this would 
involve also neglecting the depletion of the pair poten- 
tial in S. However, such assumptions, which would lead 
to the proximity effect being for many purposes inde- 
pendent on the choice of N material, could hardly be 
justified, theoretically or experimentally. 

In this work, we have carefully and rigorously exam- 
ined the various implications that the influence of a re- 
pulsive effective electron interaction in the N region has 
of the proximity effect in an N/S bilaycr. In addition to 
the spatial variation of the pair amplitude, one also has 
to study the decay of the finite pair potential in N, away 
from the N/S interface, and its depletion in S. Each 
of those spatial dependences are strongly affected both 
by the effective interaction in TV and by the N/S inter- 



facial properties. In the N material, they have oppo- 
site trends in the magnitude: F(z) is suppressed while 
|A(z)| is enhanced by a stronger repulsive interaction. 
In the superconductor, however, they are both depleted 
in the same way. This suppression of the pair potential 
near the interface is another signature of the proximity 
effect. It directly depends on the strength of this repul- 
sive interaction. We also consider the dependence of the 
proximity effect on interface scattering. Many studies 
of superconducting junctions employ a single-parameter 
description (for example, using the corresponding normal 
state reflection or transmission coefficient) for the inter- 
face properties. We show explicitly that this is clearly in- 
sufficient, even for the commonly used S- function model 10 
of the interfacial barrier. More specifically, the nature of 
the proximity effect changes independently with both the 
strength of the interfacial barrier Hb and the Fermi wave 
vector mismatch A between the N and S regions. 

We have not found any support for the formation of 
a zero bias (or near zero bias) peak in the local den- 
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sity of states (LDOS) near the N/S interface, usually at- 
tributed to the formation of Andreev bound states. We 
do find a plethora of in-gap states attributable to these 
bound states in agreement with previous work. We find 
also similar zero E features for the differential LDOS, 
after performing a subtraction of the LDOS for vanish- 
ing repulsive interaction. Such differential LDOS peaks 
become more pronounced with increasing repulsive in- 
teraction and resemble the zero bias conductance peaks 
studied extensively in unconventional superconductors. 

A challenge for future work would be to identify spe- 
cific materials and systems where the explored proximity 
effects could be readily observed. Our results show that, 
as might have been expected, the effects of a repulsive 
interaction are quantitatively important but not qualita- 
tively obvious: there is no simple and evident "smoking 
gun" . The peaks we find at small energies are in SN(E) 
not in the local DOS itself. However, careful quantita- 
tive studies of the LDOS near the interface as a func- 
tion of thickness should be (see e.g. Figs. H] through 
HU1) revealing. Such studies have been technically pos- 
sible for several years and have been found useful (see 
e.g. Refs. r50U60h in the study of Fj S systems. Further- 
more, a repulsive interaction could be considered as a 
simple model for a strongly-correlated N region. Mate- 
rials such as VO or Pd could be used as the material 
forming the N layer, as well as other materials that may 
have an enhanced susceptibility, close to the Stoner in- 
stability. Such materials should have a repulsive effec- 
tive interaction. Another direction would be to further 
examine semiconductors as the N region. Two classes 
of materials could be suitable candidates: ferromagnetic 
(III,Mn)V semiconductors, which have revealed unusual 
Andreev reflection in N/S junctions; 43 ' 66 ! 67 and nonmag- 
netic narrow-band gap semiconductors. In the second 
class, it might be useful to focus on InAs-based semicon- 



ductors. These materials offer high mobility and a sup- 
pressed Schottky barrier with S region and are already 
known for intriguing properties of Andreev reflection and 
proximity effects 6 -?. Gating of such a two-dimensional 
semiconductor would offer a natural path to alter the 
strength of the repulsive interactions. 

Finally, recent experimental and theoretical advances 
could be used in the future to extend previous ideas about 
employing screening effect o 25 ' 26 ' 28 to extract the strength 
of repulsive interactions in the N region. With an ap- 
plied magnetic field the proximity induced superconduc- 
tivity in the N region implies that there will also be a 
supercurrent and a Meissner effect. In the linear ap- 
proximation for the Meissner effect (with the supercur- 
rent proportional to the superfluid velocity) , the presence 
of 7r-states was associated with paramagnetic instabil- 
ity at N/S interface.— However, it has been shown that 
the region with suppressed pair potential can also lead 
to an important nonlinear contribution to the Meissner 
effect i 65 ' 69 Revisiting the interplay of the screening re- 
sponse on the proximity effects could provide additional 
insights for probing repulsive interactions. Thus, there is 
reason to expect that additional ways of obtaining and 
interpreting experimental data in order to extract the 
effective pairing interaction will soon be available. 
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